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The diffraction limit for lensless imaging, 
defined as the sharpest possible point im- 
age obtainable with a pinhole aperture, is 
analyzed and compared to the corre- 
sponding limit for imaging with lenses by 
means of theoretical considerations and 
numerical computations using the Fresnel- 
Lommel diffraction theory for circular 
apertures. The numerical result (m = tt) ob- 
tained for the best configuration parame- 
ter u which defines the optical setup is con- 
sistent with the quarter-wave criterion, 
and is the same as the value reported in a 
classical paper by Petzval but smaller 
than the value {u = I.Stt) found by Lord 
Rayleigh. The smallest discernible detail 
(pixel) in a composite image is defined by 
an expression found by Rayleigh on ap- 



plying the half-wave criterion and is shown 
to be consistent with the Sparrow crite- 
rion of resolution. The numerical values of 
other measures of image size are re- 
ported and compared to equivalent parame- 
ters of the Fraunhofer-Airy profile that 
governs imaging with lenses. 
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1. Introduction 



Lensless imaging is known as the basis of the camera 
obscura, or pinhole imagery. It has modern technical 
and scientific applications where the use of lenses is not 
possible (e.g., x-ray optics) or must be avoided (e.g., 
illumination of absolute radiometers). The underlying 
concept is illustrated in Fig. 1: A small circular aperture 
sd limits the pencil of rays proceeding from a luminous 
point Po in an object plane ^o and, thus, produces an 
image in the form of a small patch of light BPCB' on a 
screen 5^. It is customary to describe this image in terms 
of a dimensionless image coordinate v and configuration 
parameter u defined as 



V = - 



lirac 

IT 



U = ' 



27T a\ro + r) 
Xror 



(1) 



where A is the wavelength of light, a is the radius of the 
aperture, ro and r are the object and image distances 
shown in the figure,^ O is the aperture center, A is a 
point on the aperture rim, and c = CP is the radial dis- 
tance from the image center C. As it is the object of this 
paper to assess the actual width of the image produced 
by an aperture of given radius a, we will use v/Vu 
instead of v as the image coordinate because it is propor- 
tional to c but independent of a . 



A paraxial approximation, cosjS = 1, is assumed throughout this pa- 
per because this assumption is implicit in Eqs. (4a,b) and (6a). Without 
it, ro and r would be equal to the distances PqO and OP, respectively, 
and inclination factors would appear in most equations. 
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Fig. 1. Lensless imaging by a circular aperture M of diameter 2a. ^o = object plane, Pq = object point, 
O = aperture center, 5^= screen, BPCB' = geometrical image, c = PC = radial distance from image center. 



A general idea of the nature of the diffraction Hmit for 
lensless imaging can be obtained from a classical paper 
by Petzval [1], who assumed that for Tq » r every im- 
age point P is spread by Fraunhofer diffraction and thus 
has a finite width on the order of 2Xr/a . The aggregate 
image width 2Ac exceeds the geometrical width 2a by 
the same amount so that Ac = a -\- Xr/a, or Av/vu as 
shown by the upper curve in Fig. 2. This total image 
radius has a minimum for Vm = tt/Vw so that 

(W)petzval = TT, (Av/Vw)petzval = 2VtT, 



where (2a )petzvai is the aperture diameter expected to give 
the sharpest image. Petzval's analysis is admirable on 
account of its simplicity, although the Fraunhofer ap- 
proximation is not applicable for lensless imaging and 
the incoherent superposition of elementary diffraction 
patterns is inadmissible. 

It is easy to show that Petzval's result for u is consis- 
tent with the quarter-wave criterion of resolution. Let A 
be a point located on the spherical wavefront "^T incident 
on the aperture and on a marginal ray passing at the 
aperture center, so that 



(2(3 )p, 



ro»r, 



(2) 



AC^ = PoA^ + PoC - 2PoAPoCcosy, 



(3a) 



iS 

M^^ 




'm'- 



Fig. 2. Petzval's estimate. The diffraction limit occurs at the intersec- 
tion, « = IT, of the image widths due to diffraction (left) and geometri- 
cal optics (right). 



as may be seen by applying the cosine theorem to the 
triangle PqAC. In the paraxial approximation cosjS -- 1 
one finds PqA = PqO ^ Vq and AC -- r+ Ar, where Ar is 
the small path difference AC — OC. Therefore, 



Ar- 



^o(ro + r)(l - cosy) 



aHj^^>^ (3b) 



where it was assumed that Ar^ « r^ and y was evalu- 
ated as the small angle a fro. This proves that the quarter- 
wave criterion is satisfied when w = it. 

A further refinement of the theoretical treatment re- 
quires the Fresnel-Lommel equations [2] which govern 
the diffraction effects of lensless imaging, 



£'(m,v) =£'geom(v)la(w,v)P, £'geom(v) = 



^0^0^ 



(4a) 
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\a(u,v)P 



-[L\u,v)-\-M\u,v)l 



(4b) 



where E(u,v) is the irradiance at the point P of Fig. 1, 
£^geom(v) is the geometrical irradiance in the absence of 
diffraction, ^o is the radiant flux admitted by the aper- 
ture, \a(u,v)P is the modification of ^geom(v) by diffrac- 
tion, and L(w,v) and M(w,v) are functions defined by 
Lommel as linear combinations of infinite series of Bes- 
sel functions. 

In the past, computations based on these equations 
were tedious but nonetheless Lommel provided numeri- 
cal tables of E(u,v) and Rayleigh [3] used these tables 
for a further analysis of pinhole imaging. Rayleigh plot- 
ted the diffraction profiles which are reproduced here as 
Fig. 3 and, without additional calculations, judged that 
"w = V2'TT is too large and u = Si: is too great. The only 
question that can arise is between u = n and u = 2ii. 
The latter has decidedly the higher resolving power, but 
the advantage is to some extent paid for in the greater 



diffusion of light outside the image proper " He con- 
ducted visual and photographic experiments to settle 
this question and found that the sharpest images were 
obtained for 

WRayieigh =1.817, (2a)R^yieigh = V3.6Aror/(ro + t). (5) 

Rayleigh did not explicitly state the image size corre- 
sponding to this value of u . However, he showed that, for 
any aperture diameter 2a and in the limit ro » r, the 
greatest path difference at the point P in Fig. 1 is 2ac/r, 
so that ''the illumination will not be greatly reduced 
until the extreme discrepancy of phase reaches half a 
wavelength.'' This gives 2Ac = Xr/2a or Av = O.Stt, 
which turns out to be an excellent estimate of the onset 
of resolution in the sense of the Sparrow criterion [4] 
because the relative central irradiance^ in the com- 
posite pattern of two Lommel profiles separated by 
2Av/\4 = 'tt/\4, 





Fig. 3. Normalized Fresnel -Lommel diffraction profiles for u = 0.5 it, it, 217, and Sir. 



^ Equation (6a) is approximate in that it ignores the "shrinkage" 
caused by the fact that the maxima of the composite pattern do not 
coincide exactly with the maxima of the superimposed profiles. 
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1 +E{u,2^vl^fu) 



(6a) 



is 0.9974 for m = it, 0.9967 for u = I.Stt, and 1.0086 for 
the Airy profile defined by Eq. (7a), below. Accord- 
ingly, the smallest discernible detail (pixel) in a com- 
posite image is given by 



AV; 



Sparrow ' 



: 0.517. 



(6b) 



2. Analysis 

As the use of Lommel's equations [Eqs. (4a,b)] for 
numerical computations is no longer a problem, it is the 
purpose of this paper to re-examine the question of an 
optimal aperture diameter for lensless imaging and to 
characterize the resulting image distribution of irradi- 
ance and radiant flux in quantitative terms. This work 
was performed using the algorithms for Fresnel diffrac- 
tion published in Ref. [5], which should be consulted for 
mathematical and computational details. Readers wish- 
ing to perform computations may also consult a recent 
paper on a similar topic by Shirley [6]. 

When a lens is placed at the aperture M in Fig. 1 and 
0^0 and 5^ are conjugate object and image planes, the 
irradiance distribution in the diffraction pattern will be 
given by the Fraunhofer-Airy formula. 



E(v) = 



CpoTTfl^ 


r2j,(v)i 

V 


AV^ 



(7a) 



where Ji is the first-order Bessel function of the first 
kind and the other quantities are the same as above. As 



this equation also applies to Fresnel-Lommel diffraction 
in the limit u — > 0,^ it will be useful to analyze the 
profiles in Fig. 3 by criteria which are equivalent to 
accepted criteria for assessing lens images. The proper- 
ties of the Airy function [Eq. (7a)] that will be used for 
this purpose are listed in Table 1, where the values in the 
last column represent the fractions of the total flux ^o 
contained in a circle of radius Av , 



/(Av) 



Av 



2Ji(v) 



= 1 -J?(Av)-J?(Av), (7b) 



according to a formula derived by Rayleigh in an article 
on the wave theory of light [6] . The best known quantity 
shown in Table 1 is the radius AvAiry of the central 
maximum (Airy disk). Because 84 % of the total flux is 
concentrated in the central disk, it is often assumed that 
the diffraction pattern consists of the Airy disk alone. 
Figure 3 shows at a glance that width criteria based 
on a single irradiance value would be unreliable mea- 
sures of image sharpness as they favor the central por- 
tion of the profiles and ignore the formation of a second 
maximum which is evident in the figures. For example, 
the halfwidth of the profiles decreases monotonically 
until, near u = 317, the profile turns up before reaching 
the half-power point. On the other hand, the area width 
defined in the third row of Table 1 is quite suitable as it 
measures the feature which is most obvious on visual 
inspection of these profiles; namely, the concentration 
of E(u,v)/E(ufi) near the ordinate axis. Rayleigh may 
have used it intuitively when judging these profiles, and 
in this work it was applied mathematically by means of 
the quadrature formula 



Table 1. Properties of the Airy diffraction profile [Eq. (7a)] 



Property 


Definition 


Av 


E{Av)/E{0) 


/(Av) 


Sparrow limit 


Av = 0.5 IT 


0.5 17 


0.520 855 


0.455 925 


Area width 


Av = dv 




r2Ji(v)i 

V 


2 .^ 

:P- = 0.540 38017 

317 


0.463 082 


0.506 873 


Rayleigh width 


Av = 0.5 AvAiry 


0.609 83517 


0.367 516 


0.588 443 


Airy disk, AvAiry 


E(v) = 


1.219 67017 





0.837 785 















^ In this limit, Eq. (7a) results because one has \a(u,v)f' -^ (wJi(v)/v)^. 
However, it is misleading to refer to the Fraunhofer-Airy formula as 
a "pinhole Umit". The aperture diameters used in pinhole photography 
are on the order of 0.5 mm to 1 mm and, as shown by Eqs. (2) and 
(5), correspond to fairly large values of u except for infinitely large 
image distances. 
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1 



u\a{ufi)\ 



; I dvlQ;(w,v) 



8v 



wlaol «=o 



^\an\\ (8) 



which allows computations of area widths by straight- 
forward summation of Iq;„P = Iq;(m,v„)P for equidistant 
arguments v^ = (/i + V2) 8v up to a largest value A^ for 
which la^P is suitably small. The increments and upper 
limits used were 8v = 0.01 u and A^= 500, resulting in 
numerical values accurate to six digits or better. As 
shown in Fig. 4, these computations revealed the exis- 
tence of a shallow minimum near 



u = 1.777, Av/V^ = 0.910 925, 



(9) 



which differs insignificantly from the value I.877 found 
by Rayleigh, and thus verifies that the latter is consistent 
with the area width criterion. 



u/n 



Fig. 4. Area width Av/Vm and flux fractions /(Av/Vm) as functions 
of u. 



However, what really determines the sharpness of the 
image is not the concentration of irradiance as plotted in 
Fig. 3 but the physical concentration of radiant flux in 
the image plane itself. This aspect was mentioned but 
not explored by Rayleigh, and will be considered next. 
The quantity required for this purpose is the fractional 
flux contained in a circle of radius Ac = XrAv/lira in 
the image plane, as in Eq. (6b), and thus Eqs. (4a,b) 
were used to derive the following summation formula 
for Lommel profiles, 

Ac Av 

f(Av/Vu) = ^1 dc cE(u,v) = j^l dv v\a(u,vf 


- f^2(2^ + i){i««i' + l«-il'}. (10) 



where the notation is the same as in Eq. (8) but v and 
\a(u,v)P are replaced by their arithmetic means for each 
element of summation. Equation (10) was used to gener- 
ate lists, accurate to six digits, of/(Av/ Vw) for consec- 
utive upper limits A^ ^ 500 and increments 8v = 0.01 w. 
These lists were used as lookup tables, and linear inter- 
polation was used to obtain final results for given values 
of u and Av/Vw. 

The values of /(Av/ Vm) obtained in this manner for 
the area widths computed earlier are plotted as the lower 
curve in Fig. 4. While approaching 50 % for small val- 
ues of w , as should be expected from Table 1 , these flux 
fractions decrease for larger values of u and, where 
Av/ Vw is a minimum, they are only on the order of 
35 %. This was deemed insufficient forjudging the pro- 
files in their entirety. Therefore, it was decided to apply 
Eq. (10) in a different manner so that it would directly 
yield the values of u for which /(Av/Vw) is contained in 
the smallest possible width Av/vu). The flux fractions 
used for these computations were 0.3333, 0.5, 0.6667, 
0.75, and 0.8378, the latter being equivalent to the Airy 
disk in Table 1. The results obtained are illustrated in 
Fig. 5, showing that a shallow but discernible minimum 
of Av/vu exists for every /(Av/Vm). The values of u at 
which these minima occur are all different. However, as 
shown by Fig. 6, they are neatly clustered around a 
median value of 



w = 77, 



(11) 



which was considered the best overall choice of u based 
on the concentration of radiant flux in the image. 



Av/Vu 




Fig. 5. Flux widths Av/V^ for /(Av/VI^ = 0.3333, 0.6667, and 
0.8378. 
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25 5 75 f(Av/Vu) ^ 

Fig. 6. Dependence of best u on flux fraction /(Av/Vm). 

As a final test, we calculated the relative central irra- 
diance using Eq. (6a) in the composite pattern of two 
Lommel profiles separated by given amounts 2Av/ Vw. 
The results obtained are plotted in Fig. 7, showing that 
for the usual definition of resolving power in terms of a 
relative central irradiance on the order of 0.7 the config- 
uration parameter m= I.Stt is superior. On the other 
hand, the advantage lies with w = tt when larger separa- 
tions are considered, as has often been advocated. 



1 - 


V 






1 


\ \ ^~^ 






0.5- 


^ 


u=1.8Tt 




n- 


^~ — H \ \ 1 



0.5 
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Fig. 7. Relative central irradiance, E^JE^^^, in the composite pattern 
of two Lommel profiles separated by 2Av/Vm for m = it and I.Stt. 



3. Conclusions 

On account of the complicated nature of the Lommel 
profiles it is no surprise that the above analysis gives no 
unequivocal answer regarding a best configuration 
parameter u for lensless imaging. We have verified 
Rayleigh's value as a minimum of area width, con- 
firmed Petzval's value by computations of flux widths, 
and were unable to make a choice using the theory of 
resolution. This ambiguity can undoubtedly be at- 
tributed to the fact that the various minima of image size 



found in this analysis are all shallow so that, on the 
whole, the difference between w = tt and I.Stt is in- 
significant for practical purposes. The same can be in- 
ferred from the observation made by writers on pinhole 
photography that, although the best aperture diameters 
are usually stated within 0.01 mm, deviations on the 
order of 0. 1 mm have little effect on image quality. It 
appears that Rayleigh reached a similar conclusion. He 
mentioned that Petzval's value, w = tt, was quoted in a 
pamphlet on pinhole photography and remarked that the 
corresponding '^detail in a photograph ... was not 
markedly short of that observable by direct vision. " At 
the same time, he stated that images obtained for 
w = I.Stt 'fully bore out expectations. " 

It would of course be possible to compromise and 
adopt an intermediate best value of w. This seemed 
inadvisable as it might be misinterpreted as an improve- 
ment or refinement of the results obtained by Petzval 
and Rayleigh. The above analysis has shown that both 
are valid and that there is a fairly wide range of w's 
which give acceptably sharp images. This is in fact a 
practical advantage, as it permits the use of one and the 
same aperture for a variety of image distances and/or 
wavelengths. 

However, from a theoretical point of view it is desir- 
able to have a value which is, not only consistent with 
experience, but also conforms to accepted criteria used 
elsewhere for assessing image quality. On this basis, 
w = TT is preferable as it is consistent with the quarter- 
wave criterion whereas m= I.Stt is not. By the same 
token, the Sparrow width in Eq. (6b) is theoretically 
superior to the value in Eq. (2) because it conforms to 
the half-width criterion. Accordingly, the final result 
adopted in this paper for the diffraction limit of lensless 
imaging is 



M = TT, 



Av/V^ ^ 0.5V^, 2a = V2Aror/(ro + r), 

(12) 



where 0.5Vtt is the pixel size, defined as the smallest 
discernible detail in accordance with the Sparrow crite- 
rion. Equation (12) should not be construed as a criti- 
cism of Rayleigh's astute assessment of the Lommel 
profiles in Fig. 3. It merely ensures that the diffraction 
limits for imaging with and without lenses are defined 
consistently, but does not affect the practical aspects of 
Rayleigh's work. Yd, from a historical perspective it is 
interesting that Rayleigh did not apply his own quarter- 
wave criterion to define u while, at the same time, he 
used a half-wave criterion to define the image size. 

The optical and radiometric properties of the Lom- 
mel profile for M = TT are summarized in Table 2. Four of 
these parameters were defined to be equivalent to the 
values given in Table 1 for Airy profiles. The fifth is the 
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Table 2. Properties of the Lommel Diffraction Profile [Eqs. (4a,b)] for m = it 



Property 


Definition 


Av/V^ 


£rcl 


/(Av/\4) 


^min'^max 


Sparrow limit 


Av = 0.577 


0.886 227 


0.525 145 


0.380 803 


0.997 362 


Area width 


Eq. (7) 


1.023 941 


0.417 900 


0.458 319 


0.808 131 


Rayleigh width 


Av = 0.5 AvAiry 


1.498 744 


0.136 878 


0.647 732 


0.265 165 


Geometrical width 


AV^TT 


1.772 454 


0.060 510 


0.698 182 


0.123 190 


Airy disk, AvAiry 


/(Av/\4) = 0.837 785 


2.997 488 


0.032 497 


0.837 785 


0.064 741 



geometrical width Av/ Vm = vu that would be obtained 
in the absence of diffraction."^ The remaining columns of 
Table 2 show the numerical values of Av/ Vw and the 
corresponding relative irradiances Erc\ = E(u,Av/vu)/ 
E{ufi), flux fractions /(Av/Vm), and relative central 
irradiance E^JEmax in a double image. The equivalent 
Airy disk in Table 2 is defined so that the corresponding 
flux fraction is the same as for the Airy profile, and the 
Rayleigh width is half as large;^ their numerical values 
can be closely approximated by 3 and 1.5, respectively. 
It should be noted that the geometrical width is twice as 
large as the Sparrow width and appears near the bottom 
of the table. This shows that, on the whole, the diffrac- 
tion-limited image is sharper than the geometrical im- 
age. 
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^ This quantity is germane to Lommel profiles and is zero in the Airy 
limit M ^ 0. 

^ The relative central irradiance £min/£^max corresponding to this 
Rayleigh width is only 0.27. Alternatively the Rayleigh width of the 
Lommel profile can be evaluated as Av/Vm= 1.075 287, which is 
smaller than the value in Table 2 but gives the familiar relative central 
irradiance 0.735 932 found for Airy profiles. 
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